

**********************************************************************************
*** Spatial Inequalities in Premature Mortality in Great Britain 
*** Authors in alphabetical order: Denise Laroze, Eric Neumayer, Thomas Pluemper.
*** STATA 12 
*** Replication data and additional non-reported results
********************************************************************************** 



***********************************************************************************
local DIR = "D:\Research\Health for all database\WHO mortality database\"
  /*change relative path to the directory where the files are located */
cd "`DIR'"
***********************************************************************************




use "ipmgb.dta", replace

set more off



***********************************
******** Detecting Spatial Patterns
***********************************



********************
******** Table 1
********************


** Great Britain
* Male 
reg pm0_70_m    northerness westerness, robust
outreg2 using table1, excel 2aster  replace  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f    northerness westerness, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg pm0_70_m    northerness westerness if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f    northerness westerness if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg pm0_70_m    northerness westerness if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f    northerness westerness if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 



** Great Britain
* Male 
reg pm0_70_m    centrality, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg pm0_70_f    centrality, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg pm0_70_m    centrality if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f    centrality if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg pm0_70_m    centrality if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f    centrality if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 



** Great Britain
* Male 
reg pm0_70_m     pm070_m_splag_cont_rwst, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg pm0_70_f     pm070_f_splag_cont_rwst, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg pm0_70_m     pm070_m_splag_cont_rwst if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f     pm070_f_splag_cont_rwst if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg pm0_70_m     pm070_m_splag_cont_rwst if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f     pm070_f_splag_cont_rwst if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 




** Great Britain
* Male 
reg pm0_70_m     pm070_m_splag_prox_rwst, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg pm0_70_f     pm070_f_splag_prox_rwst, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg pm0_70_m     pm070_m_splag_prox_rwst if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4) label ctitle(PM Male) 

* Female 
reg pm0_70_f     pm070_f_splag_prox_rwst if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg pm0_70_m     pm070_m_splag_prox_rwst if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f     pm070_f_splag_prox_rwst if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 




** Great Britain
* Male 
reg pm0_70_m    area_density , robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg pm0_70_f    area_density, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg pm0_70_m    area_density if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f    area_density if country==1, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg pm0_70_m    area_density if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f    area_density if country==3, robust
outreg2 using table1, excel 2aster  append  dec(4)  label ctitle(PM Female) 




***********************************
******** The Socio-Economic Determinants
***********************************



********************
******** Table 2
********************

preserve

* Express income as income in �k
replace mean_income_imp=mean_income_imp/1000

* Express premature mortality as probability
replace pm0_70_m=pm0_70_m/100000

* Express premature mortality as probability
replace pm0_70_f=pm0_70_f/100000

* Descriptive summary stats of premature mortality
su pm0_70_m
su pm0_70_f

*** Male
reg pm0_70_m mean_income_imp benefit_claimants_prop_1664 /*
	 */ prop_q_L1 prop_q_L2  prop_q_L3 prop_q_L4 /* 
	 */prop_ind_agri prop_ind_mining prop_ind_manufacture prop_ind_elect_gas prop_ind_water prop_ind_construction /*
	 */prop_ind_retail prop_ind_transport prop_ind_hotel_food prop_ind_information prop_ind_financial prop_ind_real_estate /*
	 */prop_ind_profess_scientific prop_ind_admin prop_ind_public_admin prop_ind_education prop_ind_health /*
	 */ prop_L1_sec prop_L2_sec prop_L3_sec prop_L4_sec prop_L5_sec prop_L6_sec prop_L7_sec/*
	 */ prop_mixed prop_asian prop_black prop_white,  
outreg2 using table2, excel  replace noaster  dec(4)  label ctitle(PM Male)


*** Female
reg pm0_70_f mean_income_imp prop_L1_sec prop_L2_sec prop_L3_sec prop_L4_sec prop_L5_sec prop_L6_sec prop_L7_sec/*
	 */ prop_mixed prop_asian prop_black prop_white  /*
	 */ prop_q_L1 prop_q_L2  prop_q_L3 prop_q_L4 /* 
	 */ benefit_claimants_prop_1664 /*	
	 */prop_ind_agri prop_ind_mining prop_ind_manufacture prop_ind_elect_gas prop_ind_water prop_ind_construction /*
	 */prop_ind_retail prop_ind_transport prop_ind_hotel_food prop_ind_information prop_ind_financial prop_ind_real_estate /*
	 */prop_ind_profess_scientific prop_ind_admin prop_ind_public_admin prop_ind_education prop_ind_health,  
outreg2 using table2, excel  append noaster dec(4)  label ctitle(PM Female)

restore




***********************************
******** Explaining the Spatial Patterns
***********************************



*** Male
quietly reg pm0_70_m mean_income_imp benefit_claimants_prop_1664 /*
	 */ prop_q_L1 prop_q_L2  prop_q_L3 prop_q_L4 /* 
	 */prop_ind_agri prop_ind_mining prop_ind_manufacture prop_ind_elect_gas prop_ind_water prop_ind_construction /*
	 */prop_ind_retail prop_ind_transport prop_ind_hotel_food prop_ind_information prop_ind_financial prop_ind_real_estate /*
	 */prop_ind_profess_scientific prop_ind_admin prop_ind_public_admin prop_ind_education prop_ind_health /*
	 */ prop_L1_sec prop_L2_sec prop_L3_sec prop_L4_sec prop_L5_sec prop_L6_sec prop_L7_sec/*
	 */ prop_mixed prop_asian prop_black prop_white, robust 
predict resi_reg1_m, resid
predict y_star0_70_m // used for the maps
predict resid0_70_m, resid // used for the maps


*** Female
quietly reg pm0_70_f mean_income_imp prop_L1_sec prop_L2_sec prop_L3_sec prop_L4_sec prop_L5_sec prop_L6_sec prop_L7_sec/*
	 */ prop_mixed prop_asian prop_black prop_white  /*
	 */ prop_q_L1 prop_q_L2  prop_q_L3 prop_q_L4 /* 
	 */ benefit_claimants_prop_1664 /*	
	 */prop_ind_agri prop_ind_mining prop_ind_manufacture prop_ind_elect_gas prop_ind_water prop_ind_construction /*
	 */prop_ind_retail prop_ind_transport prop_ind_hotel_food prop_ind_information prop_ind_financial prop_ind_real_estate /*
	 */prop_ind_profess_scientific prop_ind_admin prop_ind_public_admin prop_ind_education prop_ind_health, robust 
predict resi_reg1_f, resid
predict y_star0_70_f // used for the maps
predict resid0_70_f, resid // used for the maps





********************
******** Table 3
********************


** Great Britain
* Male 
reg resi_reg1_m    northerness westerness, robust
outreg2 using table3, excel 2aster  replace  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f    northerness westerness, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg resi_reg1_m    northerness westerness if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f    northerness westerness if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg resi_reg1_m    northerness westerness if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f    northerness westerness if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 



** Great Britain
* Male 
reg resi_reg1_m    centrality, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg resi_reg1_f    centrality, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg resi_reg1_m    centrality if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f    centrality if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg resi_reg1_m    centrality if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f    centrality if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 



** Great Britain
* Male 
reg resi_reg1_m     pm070_m_splag_cont_rwst, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg resi_reg1_f     pm070_f_splag_cont_rwst, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg resi_reg1_m     pm070_m_splag_cont_rwst if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f     pm070_f_splag_cont_rwst if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg resi_reg1_m     pm070_m_splag_cont_rwst if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f     pm070_f_splag_cont_rwst if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 




** Great Britain
* Male 
reg resi_reg1_m     pm070_m_splag_prox_rwst, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg resi_reg1_f     pm070_f_splag_prox_rwst, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg resi_reg1_m     pm070_m_splag_prox_rwst if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4) label ctitle(PM Male) 

* Female 
reg resi_reg1_f     pm070_f_splag_prox_rwst if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg resi_reg1_m     pm070_m_splag_prox_rwst if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f     pm070_f_splag_prox_rwst if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 




** Great Britain
* Male 
reg resi_reg1_m    area_density , robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg resi_reg1_f    area_density, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** England
* Male 
reg resi_reg1_m    area_density if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f    area_density if country==1, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 


** Scotland
* Male 
reg resi_reg1_m    area_density if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f    area_density if country==3, robust
outreg2 using table3, excel 2aster  append  dec(4)  label ctitle(PM Female) 





********************
******** Table 4
********************

capture dropvars westerness
gen westerness=-easterness


** Great Britain
* Male 
reg pm0_70_m    northerness westerness, robust
outreg2 using table4, excel 2aster  replace  dec(4)  label ctitle(PM Male) 

* Female 
reg pm0_70_f    northerness westerness, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Male 
reg resi_reg1_m    northerness westerness, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Male) 

* Female 
reg resi_reg1_f    northerness westerness, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 







** Great Britain
* Male 
reg pm0_70_m    centrality, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg pm0_70_f    centrality, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 


* Male 
reg resi_reg1_m    centrality, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg resi_reg1_f    centrality, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 




** Great Britain
* Male 
reg pm0_70_m     pm070_m_splag_cont_rwst, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg pm0_70_f     pm070_f_splag_cont_rwst, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 


* Male 
reg resi_reg1_m    pm070_m_splag_cont_rwst, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg resi_reg1_f    pm070_f_splag_cont_rwst, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 





** Great Britain
* Male 
reg pm0_70_m     pm070_m_splag_prox_rwst, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg pm0_70_f     pm070_f_splag_prox_rwst, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 


* Male 
reg resi_reg1_m    pm070_m_splag_prox_rwst, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg resi_reg1_f    pm070_f_splag_prox_rwst, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 





** Great Britain
* Male 
reg pm0_70_m    area_density , robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg pm0_70_f    area_density, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 


* Male 
reg resi_reg1_m    area_density, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 

* Female 
reg resi_reg1_f    area_density, robust
outreg2 using table4, excel 2aster  append  dec(4)  label ctitle(PM Female) 




********************
******** The Glasgow Effect
********************

** Male excess mortality in Glasgow
su pm0_70_m
di r(max)-r(mean)

** Female excess mortality in Glasgow
su pm0_70_f
di r(max)-r(mean)

** Male excess UNEXPLAINED mortality in Glasgow
preserve
keep if area==157
di resi_reg1_m 
restore

** Female excess UNEXPLAINED mortality in Glasgow
preserve
keep if area==157
di resi_reg1_f 
restore


